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Abstract 

There is a well established theory of Ballooning modes in a toroidal plasma. The 
cornerstone of this is a local eigenvalue A on each magnetic surface - which also depends on 
the Ballooning phase angle k. In stationary plasmas, X(k) is required only near its 
maximum, but in rotating plasmas its average over k is required. Unfortunately in many 
cases X(k) does not exist for some range of k, because the spectrum there contains only a 
stable continuum. This limits the application of the theory, and raises the important 
question of whether this "stable interval" gives rise to significant damping. This question is 
re-examined using a new, simplified, model - which leads to the conclusion that there is no 
appreciable damping at small shear flow. In particular, therefore, a small shear flow should 
not affect Ballooning mode stability boundaries. 



1. Introduction 

MHD Ballooning modes, that is perturbations with short wavelength perpendicular to the 
local magnetic field, but long wavelength parallel to it, are a particularly insidious form of 
plasma instability pp. In an axi-symmetric toroidal plasma they are normal modes with 
large toroidal mode number n. If the plasma is stationary (no toroidal flow) then after 
introducing an extended poloidal coordinate rj (which removes explicit periodicity 
constraints [2]) these modes can be represented by an eikonal that describes the short wave 
behavior, 

£ = rj) expinq'(xr) + S(x)) (1) 

where x is a minor radius coordinate, constant on each magnetic surface, and q'(x) is the 
average magnetic shear. 

In the limit n — > oo, each magnetic surface then decouples from its neighbors and has its 
own local growth rate X(x, k), given as the eigenvalue of an ODE in rj, with the surface label 
x and the phase k = dS/dx (related to the poloidal location of the mode) as parameters [2], 

^P(v, x, k) — 77) + Q(rj; x, k) 77) = X 2 (x, k) R(r]; x, k) £(77) (2) 

The coefficients P,Q,R, are defined by the plasma equilibrium and are periodic in the 
phase angle k. The boundary condition is that £ be bounded as \q\ — >■ 00 and the 
eigenvalue X(x, k) is also periodic in k. 

This local growth rate X(x, k) contains the information necessary to construct global 
Ballooning modes [2]. For example, near marginal stability, the Ballooning mode growth 
rate is 
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where X max is the maximum of X(x, k) wrt both x and k. 



If there is a sheared toroidal rotation Q(x) in the equilibrium plasma, the eikonal (1) is no 
longer valid. Instead the appropriate form is [3] 

£(x,r)) exp[inq'(xr] + S(x) + fl(x)t)} (4) 



In the limit n — > oo, each magnetic surface again decouples from its neighbors but its time 
dependence is no longer purely exponential. It is of Floquet form, given [3] by the PDE 

~P( V ; x, k) ~ fa t) + Q(v, x, k) fa t) = ®- t R( V ; x, k) ~ fa, t) (5) 

where k = k + s v t and s v = ft'/ q' is the velocity shear. 

When the shear velocity is small the coefficients P, Q, R evolve slowly in time and we can 
find an "adiabatic" solution of EqtjiJ with exponential growth rate ~ exp(/ X(x, t)dt), 
where 

^P( V ; x, k)~ fa t) + Q( V ; x, k) fa t) = X 2 (x } k) R( m x, k) fa (6) 



Thus X(x,t) is just the local growth rate of the static plasma (given by Eq([2j)), at a time 
dependent phase angle (k + s v t), and the average growth rate of Ballooning modes in a 
rotating plasma is therefore 

A~^-{x(k)dk (7) 

So we see that, even with shear flow, Ballooning instabilities can be described through the 
local growth rate X(x, k) of a static plasma. However there is a crucial difference in the role 
of A when shear flow is significant. In the absence of shear flow, A is required only in the 
vicinity of its maximum wrt k, say at k m , whereas with shear flow it is required over the 
full range < k < 2tc. This seemingly small difference introduces a fundamental difficulty: 
X(x, k) may not be defined for all k. 



This problem arises when, as is often the case, Eq([2j) has an unstable eigenvalue only over 
a restricted range of phase angle k, around k m : outside this range it has no discrete 
eigenvalues - only a stable continuum A = ± iu>. Within this "stable interval" of k, the 
solution cannot be represented by an adiabatic form with exponential time dependence. 
(This depends on the highest eigenvalue being well separated from all others - a condition 
that is clearly violated when there is only the continuum.) It is then not clear how A 
should be interpreted (or replaced). 

Formally therefore, the applicability of Ballooning theory to plasmas with shear flow is 
severely limited. This limitation is particularly apparent in determining Ballooning mode 
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stability boundaries. By definition these are where a mode with some particular phase k 
first becomes unstable; however in general, modes with any other phase angle will remain 
stable - so that near a stability boundary we expect ballooning modes to have a significant 
stable interval. An important issue therefore is whether, with sheared rotation, this stable 
interval leads to significant damping. If not, then although even a small rotation will 
reduce the growth rate of Ballooning modes, it alone cannot stabilize them - consequently 
the formal stability boundaries for Ballooning modes will be unchanged by small rotation. 

The aim of this paper is to investigate the effect of a stable interval on Ballooning modes in 
a plasma with small velocity shear. The next Section describes a model of a toroidal 
plasma which incorporates a stable interval and introduces previous (contradictory) 
attempts to deal with it. Section 3 describes a new approach to the problem. This is 
discussed in detail in Section 4. A summary and some conclusions are given in Section 5. 

2. The stable interval 

There are conflicting views, leading to contradictory conclusions, on how to deal with the 
the stable interval. To illustrate the problem we invoke the well known (s,a) model 0] of a 
large aspect ratio Tokamak. (a is the plasma pressure gradient and s = rq'/q is the 
magnetic shear.) In this model, Eq(|5]) takes the specific form 



^ (i + p 2 ) ^ av, t)+q to, o = ^(i+ p 2 ) | av, t) (8) 

with 

P = si] — asin(?y + s v t) Q = a[cos(?7 + s v t) + Psm(r] + s v t)] (9) 
When the shear velocity s v is small, Eq(|8]) can be reduced to the simpler "wave Eq" form 

g + ^,^ = § (10) 



where ift = £/ yl + P 2 and the "potential" V i 



is 



[s - acos(?7 + s v t)} 2 a cos(r] + s v t) 
V(v,s v t) = + (i + p2) (11) 



The corresponding local Eq is 



drf 



+ V( V ,s v t)iP = \ 2 iP (12) 



A stability diagram for the s, a model is shown in Fig(l). As expected, there is a large 
region within the conventional "unstable" region where most phase angles are still stable, 
i.e where there is a significant stable interval in k. 



One way to deal with the stable interval [5] is to represent the perturbation as an integral 
over the continuum modes. There are difficulties in justifying this approach, but the 
underlying picture is that continuum modes are excited as the phase k + s v t enters the 
stable interval. The "depth" of this excitation is proportional to s v but the duration of the 
subsequent stable interval is proportional to l/s v - so that phase mixing during the stable 
interval reduces the amplitude by a fixed amount per transit through the interval. 
Averaged over the whole interval this is equivalent to a damping rate that is proportional 
to Q' and vanishes as Q' — > 0. 



Another approach [7] is based on the fact that as i] — > oo the solution of Eq(10) will consist 
of outgoing and incoming, reflected, waves. If the reflected waves can be neglected it 
becomes appropriate to impose an "outgoing wave only" boundary condition. This has no 



effect on eigenvalues of Eq(10) in the unstable interval - because if A > solutions that 
match to outgoing waves are those which also decay as rj — > oo. However, in a stable 
interval the change in boundary condition has a dramatic effect - because if A < a 
solution that matches to an outgoing wave would otherwise diverge as i] — > oo, and would 
be rejected as an eigenf unction. In fact, with the "outgoing wave" boundary condition, 



Eq(10) no longer has a continuum of eigenvalues A = ±ico in the stable interval; instead it 
will usually have a discrete negative eigenvalue. In this case the stable interval will 
introduce a significant "wave" damping that is independent of Q' and persists as O' — > . 

This is essentially the view taken by Waelbroek and Chen [7j, although they also 



introduced a further approximation. In essence they assumed that in Eq(10) the range of rj 
can be separated into an inner region where the inertial term A 2 is negligible, and an outer 



region where V is negligible. One then solves Eq(10) for i/j in the inner region (with 
A 2 = 0) and at the inner region boundary matches its logarithmic derivative A' to an 



outgoing wave. Then A = A' and the problem of finding an eigenvalue of the full Eq(10) 
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with the outgoing wave boundary condition is avoided. 



Whether or not one adopts this additional approximation, (which may be valid only for 
small A'), the crucial question is whether reflected waves can be neglected. The fact that 
the potential V is small at large rj certainly makes the reflection coefficient small. But if 
the perturbation were decaying exponentially during the stable interval, as it does when 
reflection is ignored, then any reflected wave would have been created when the outgoing 
wave was exponentially larger than it is when the reflected wave returns. This large 
exponential factor could outweigh the small reflection coefficient when s v is small - making 
the reflected wave important during the stable interval. 

3. A 'Toy' model 

To investigate the importance of reflection we introduce a simple "Toy" model. Suppose 
first that the potential vanishes outside some small region around rj — 0, that is 
V(r],t) = —D(s v t)5(r]) (so that D corresponds to A' in the account above). This Toy 
model has an exact solution 

vp = Aexp J X 1 D(s v t')dt' (13) 



valid whether D is positive (unstable interval) or negative (stable interval) and for all s v . 
This is precisely the result one would get from the Waelbroek and Chen approximation and 
appears to support their picture of the damping. But of course the model does not yet 
address the question of reflections from an extended potential. To do this we modify it by 
adding a weak potential "tail" v{rf) that — > as rj — > oo. Then V(rj,t) = D(s v t)5(i]) + v (if). 



At this point it is convenient to express Eq(10) (using the Green's function or more simply 



by introducing (rj ± t) as coordinates) in an integral form: 



iP( V , t) = \f dt' r +itt } V>(i/, V(rf, drf (14) 

2 J Jr)-(t-t') 



Then, with V(r],t) = D(s v t)8(rj) +v(rj), the central perturbation, ?/>(0,t) = ^(t) satisfies 

r/ * (/) -D(st)*(t) + l r ^rf,t-rf)v[rf,t-rf) (15) 



dt v ' v ' 2 Jo 
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Reflected waves arise from the integral term in Eq( 15). To calculate the first order reflected 
wave, linear in v, we can replace ipiv'i t ~ v') m this integral by its form when v is ignored - 
then ip(j],t) is constant along lines of fixed (t — if) and i/j(r)',t — rj') = ty(t — 2r/). Thus we 
obtain a closed Eq for the central perturbation ^(t). 

„, 1 /•* 



dt 



D(st)V(t) + - I *(v>((t ~ V')f2, (t + tj)/2)dt/ (16) 

Z J—oo 



(For a geometrical interpretation of Eqs(14 16) see Figs(2,3)). 
4. Calculation 

The effect of a weak, but extended, potential tail can be calculated explicitly in the case 
that it decays exponentially ~ exp(— prj) and is independent of time. Then 



dt 

With t = s v t, this is equivalent to the Eqs 



D(s v t)V(t) +v [ Vtf) exp(-p (t - n'))dr,' (17) 

J -co 



d^ d$ 
s v — - D{t)^(t) = $(r) : s v — + p $(r) = v ^{t) 

CLT CLT 



When s v is small, Eqs(18) have the WKB solutions 



V p (r) = q{r) exp(I f D(r')dr'), V s \r) = -L 1 exp(-^r) (19) 
where 

q(r) = exp( (^ ) [ — ) dr') (20) 
s v J {D(t') +p) 

which, since vq is assumed small, is slowly varying compared to tp p (t) , tp s (t) . 

Thus we see that perturbation \1/ has two components. One is essentially determined by the 
strength of the central potential D; the other by rate of decay (not the strength!), of the 
potential tail. For convenience we will refer to these as the Primary (P) and Secondary (S) 
components respectively. The two components are independent of each other except at the 
transition points where D + p = 0, - that is near the beginning and end of a stable interval. 

To describe the coupling between P and S components at the transition points it is 
convenient to write \P(t) = xi T ) exp(— pr/s v ) . Then 

*v^T - (D(r)+ P ) d f - (D'(r) + ^) X = (21) 
dr dr s v 



At the first transition, near the start of a stable interval, (D + p) ~ — ar with a > and 

?ir+T^r + (i-&)x = o (22) 

where u = yja/s v r, b — vo/as v and the Primary and Secondary components are 
x p {u) = u- b exp(-u 2 /2), x S (u) = u^-V. 



The general solution of Eq(22) can be expressed in terms of Parabolic Cylinder Functions as 



X = [A B_ b (u) + B D (6 _!) (w)] exp(-u 2 /4) (23) 

and the asymptotic expansions of the PCFs link the amplitudes of x P ( u ) an d X \ u ) before 
and after the transition. In particular, D __{,(«), links x F \ which is sub-dominant before the 
transition, to the combination (x P + C\ x S ) after the transition. That is 

X P (u) — ► D_ 6 (u) exp(-n 2 /4) — ► X p (u) + C lX 5 (u) (24) 

with the coupling coefficient C\ = y/2n/T(b). 



After the transition the primary component decays rapidly, so this transition effectively 
converts the exponentially growing perturbation in the unstable interval into a slowly 
varying perturbation during the stable interval. In terms of r; 

r- b exp(-ar 2 /2^) — ► d ■ {a/ s v f- 1/2) r (6_1) (25) 

At the second transition, near the end of the stable interval, D + p ~ +ar and 

**-„g-(l + ») x = (26) 
The Primary and Secondary components at this transition are x P { u ) = uh exp(w 2 /2) , 



X s = u ( b+1 ), and the general solution of (26) is 



X = [A B b (iu) + B D_ (6+1) (-«)] exp(W/4) (27) 

In this case the asymptotic expansions show that D_(& +1 )(w) links x S > sub-dominant before 
the transition, to the combination (x 5 + Ci x P ) after the transition. 

X S( U ) _^ D_ (6+1) (u) exp( M 2 /4) — ► ( X 5 ( M ) + C 2X P («)) (28) 
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with coupling coefficient C 2 = y/2n/T(l + b). After this transition the secondary 
component rapidly becomes negligible compared to the exponentially growing primary 
component. This transition therefore transforms the slowly varying perturbation during 
the stable interval back to exponential growth in the next unstable interval. In terms of r 

T -(b+i) _^ C2 . ( a /<g(Hi/2) T b exp ( ar 2 /2 ^) (29) 



To calculate the full effect of passing through a stable interval we need also to determine 
the change in the secondary component across the interval, that is from x S ( T i) = T< i ^ a t 
a time T\ immediately after the first transition, to x{ T 2) = t 2 at a time t 2 immediately 



before the second transition. From Eqs(19, 20) this is given by 



S7 \ -(6+1) 

Here fi is a numerical coefficient that depends on the precise form of D{r) in the stable 
interval. (When D(r) = — asin(r), /n = 2). 

We see therefore, that if we follow the perturbation from an unstable interval, through the 
following stable interval and into the next unstable interval, its amplitude is changed by a 
factor. C = dC 2 {fi) 2b {a/s v ) 2b . 

5. Summary and Conclusion 

There is a well established theory of MHD Ballooning modes in a toroidal plasma [U [2] . 
The cornerstone of this theory is a local growth rate X(x, k) - defined as the largest 
eigenvalue of an ODE 

_± P ( m X) k) _± + Q ( m X) k) = X \x, k) R( m x, k) C(v) (31) 

on each magnetic surface x. ( rj is an extended poloidal coordinate and k a poloidal phase; 
the boundary condition is that £(77) be bounded as 77 — > 00.) 

If the plasma is at rest X(x, k) is required only near its maximum wrt k, but if the plasma 
has a toroidal flow, the phase k increases continuously with time and A is required over the 
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full range < k < 2tt. This leads to a serious problem: in many cases A is not defined for 



all values of k. Instead there is a "stable interval" of k where Eq(31) has no discrete 
eigenvalue - only a continuum of stable eigenvalues A = iu. The existence of this stable 
interval restricts the application of Ideal MHD Ballooning theory in rotating plasmas and 
raises the question of whether it leads to significant damping - which is particularly 
important when considering stability boundaries. 

One suggestion [5] for dealing with a stable interval was to expand the perturbation in the 
(singular) continuum modes. This led to the conclusion that the stable interval produces 
only a small damping that vanishes as Q' — > 0, but it is difficult to justify this analysis. 

A seemingly more satisfactory approach [7\ is based on the fact that at large \rj\ solutions 



of Eq(31) consist of independent outgoing and incoming (reflected) waves. If the reflected 



waves can be neglected, an "outgoing wave only" boundary condition becomes appropriate 



for Eq(31). Then, in the stable interval, Eq(31) no longer has a continuum of stable 



eigenvalues: instead it may have a discrete negative eigenvalue. In this event there would 
be significant damping during the stable interval - which would be independent of Q' and 
persist as Q' — > 0. 

The crucial question in this picture is whether reflected waves can indeed be neglected. 
The analysis in this paper, based on a model that specifically considers reflections, suggests 
otherwise. It appears that no matter how small, reflections determine the behavior during 



the stable interval and the negative eigenvalues of Eq(31) are therefore not relevant. Any 
damping that does occur during the stable interval is proportional to the exponent in the 
decay of the reflection coefficient at large |?7|. 

In the (s, a) model of a Tokamak, the reflection coefficient decays only algebraically, 
implying negligible damping during the stable interval. Instead its effect is felt mainly at 
the start and finish of that interval. At the start there is a transition from exponential 
growth of the perturbation to a near constant amplitude throughout the stable interval, 
and at the end of the interval there is a transition back to exponential growth. These 
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transitions, and the stable interval, change the perturbation amplitude by a factor 
C = (27r/r(6)r(l + b)) [i 2b (a/s v ) 2b each time the phase k(t) passes through the stable 
interval. (Here s v /a determines the shear flow, v /a 2 determines the strength of reflections 
and b = vo/as v .) At small shear (b > 1) this change is equivalent to an average time 
constant throughout the interval 

A^' e = sJogC ~ ^° ( log ( — ) + 1) - (32) 
a v o v 

This result indicates that, as the shear flow Q' — > 0, [11J this damping is negligible 
compared to the growth rate in an unstable interval. It is also negligible compared to the 
damping which would occur if reflections were neglected. (In effect, therefore, the effective 
growth rate for Ballooning modes can be summarized as "A(fc) in an unstable interval and 
zero in a stable interval".) 



This behavior is shown clearly in numerical solutions of the underlying model Eqs(18). 
Figs(4,5) show that the instantaneous growth rate closely follows (exp J(D(r)/s v ) dr) 
during the unstable interval and becomes briefly negative at the start of the stable interval 
(where the coupling coefficient is less than unity). It is then almost constant during the 
remainder of the stable interval until it returns smoothly to exponential growth at the end 
of the interval (where the coupling coefficient C<i is large). Figs(6,7) confirm that as s v 
decreases the perturbation becomes effectively constant throughout the stable interval. 

These features are remarkably similar to those seen in numerical computations of 
Ballooning modes by Furukawa and colleagues 0, El [TO] . These computations simulate the 
behavior of Ballooning modes in realistic Tokamak configurations, using initial value codes. 
The example in Fig(8) (adapted from Ref.[10j), shows the logarithmic growth rate of the 
kinetic energy of an n — > oo Ballooning mode in a Tokamak with toroidal rotation. (The 
"Mach" numbers Mj are proportional to the shear flow velocity.) This closely resembles the 
model solutions in Figs(5, 7). 

In conclusion, it appears that, contrary to some previous conclusions, a stable interval in k 
does not lead to significant damping of Ballooning modes in a plasma with small toroidal 
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shear flow, and therefore does not change formal Ballooning mode stability boundaries. 

I am grateful to J Connor for many valuable discussions and comments, to J Hastie for the 
calculation of Fig(l) and to M Furukawa and N Aiba for Fig(8). 
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Fig 1 Stability Diagram for s-alpha model 
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Fig 2 Region of integration for Eq 14 
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Fig 3 Region of integration for Eqs 15,16 
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Kinetic energy (K) of infinite n ballooning mode. 
Fig adapted (with permission) from ref [10] 
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